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Abstract. We use a variational principle to prove an existence and unique- 
ness theorem for planar weighted Delaunay triangulations (with non-intersect- 
ing site-circles) with prescribed combinatorial type and circle intersection an- 
gles. Such weighted Delaunay triangulations may be interpreted as images of 
hyperbolic polyhedra with one vertex on and the remaining vertices beyond the 
infinite boundary of hyperbolic space. Thus the main theorem states necessary 
and sufficient conditions for the existence and uniqueness of such polyhedra 
with prescribed combinatorial type and dihedral angles. More generally, we 
consider weighted Delaunay triangulations in piecewise flat surfaces, allowing 
cone singularities with prescribed cone angles in the vertices. The material 
presented here extends work by Rivin on Delaunay triangulations and ideal 
polyhedra. 



1. Introduction 

1.1. Overview. Rivin developed a variational method to prove the existence and 
uniqueness of ideal hyperbolic polyhedra with prescribed combinatorial type and 
dihedral angles, or equivalently, of planar Delaunay triangulations with prescribed 
combinatorial type and circumcircle intersection angles [26]. The purpose of this 
article is to extend this method to hyperideal polyhedra and to weighted Delaunay 
triangulations. Consider a finite set of disjoint circular disks in the plane. For 
every triple of such disks there exists a circle that intersects the boundaries of 
the disks orthogonally. The weighted Delaunay triangulation induced by the disks 
consists of the triangles whose vertices are the centers of a triple of disks such that 
the orthogonal circle of this triple intersects no other disk more than orthogonally 
(see Figure 2). Such weighted Delaunay triangulations correspond to hyperbolic 
polyhedra with one vertex on the infinite boundary and all other vertices outside. 
The dihedral angles correspond to the intersection angles of the orthogonal circles. 
More generally, we allow cone singularities at the centers of the disks, and we call 
such weighted Delaunay triangulation with cone singularities euclidean hyperideal 
circle patterns. The question we consider is: Given an abstract triangulation T 
with vertex set V, edge set E, and face set T, and given intersection angles as a 
function on the edges and cone angles S as a function on the vertices, does there 
exist a corresponding euclidean hyperideal circle pattern, and is it unique? The 
main result is the following. 

Theorem 1. A euclidean hyperideal circle pattern with triangulation T , intersec- 
tion angles 6 : E [0,7i') and cone/boundary angles "E. : V —>■ (0, oo) exists if and 
only if the set of coherent angle systems A{T, 6, S) is not empty. In this case, the 
circle pattern is unique up to scale. 

Section 2 contains the basic definitions and a precise statement of the "circle 
pattern problem" under consideration. The set of coherent angle systems ■4(T, 9, S) 
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Figure 1. A Delaunay triangulation. 

is defined in Seetion 2.5. It is a subset of R^'-^l defined by a system of linear 
equations and linear strict inequalities. The proof of Theorem 1 is based on the 
variational principle that is presented in Section 3, where a function F : R^'-^l — > II 
is defined explicitely in terms of Milnor's Lobachevsky function. The critical points 
of F in A{T, 9, S) correspond to solutions of the circle pattern problem (Lemma 1). 
The uniqueness of a solution follows immediately from the fact that the function F 
is strictly concave on A{T, 9, S) (Lemma 2). To prove the existence of a solution, 
we show that F cannot attain its maximum on the boundary A^T, 9, 'E.)\A{T, 9, S) 
(Lemma 3). Sections 4-6 are devoted to the proofs of these three main lemmas. 
The explicit formula for the function F is based on a hyperbolic volume formula 
which is derived in Section 7. 

The variational principle presented here is not only a tool to prove the existence 
and uniqueness Theorem 1. Since it reduces the circle pattern problem to a convex 
optimization problem with linear constraints, it provides a means for its numerical 
solution. This is important in view of possible applications, such as using circle 
patterns to map 3D triangle meshes to the plane. 

Different necessary and sufficient conditions for the existence of hyperideal circle 
patterns were obtained by Bao & Bonahon [4] and by Schlenker [32]. We discuss 
these results in Section 1.3. 

1.2. Delaunay triangulations and hyperbolic polyhedra. "Patterns of cir- 
cles" have become objects of mathematical interest after Thurston introduced them 
as elementary and intuitive images of polyhedra in hyperbolic 3-space [35]. Thus, 
a planar Delaunay triangulation {i.e. a triangulation of a convex polygonal region 
with the property that the circumcircle of each triangle does not contain any vertices 
in its interior, see Figure 1) can be viewed as representation of a convex hyperbolic 
polyhedron with all vertices in the infinite boundary of hyperbolic space: Erase all 
interior edges, keep only the circumcircles and the boundary edges, and extend the 
boundary edges to straight lines. Consider the paper plane as the infinite boundary 
of hyperbolic space, represented in the Poincare half-space model. In this model, 
hyperbolic planes are represented by hemispheres and half-planes which intersect 
the infinite boundary orthogonally in circles and lines, respectively. Therefore, if 
we erect hemispheres and orthogonal half-planes over the circumcircles and the pro- 
longed boundary edges, we obtain a set of hyperbolic planes which bound a convex 
polyhedron. This polyhedron's vertices are the vertices of the Delaunay triangula- 
tion and one additional point, the infinite point of the boundary plane, where all 
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Figure 2. A weighted Delaunay triangulation. Note that the intersection points 
of two adjacent face-circles He on the edge between them. 



the hyperbohc planes corresponding to the boundary edges intersect. Moreover, 
the dihedral angle at an edge of the polyhedron is equal to the angle in which the 
corresponding circles/lines intersect. 

A point in the infinite boundary of hyperbolic space is called an ideal point, and 
a polyhedron with all vertices in the ideal boundary is called an ideal polyhedron. 
This terminology has become widely accepted; in the old literature, the term "ideal 
point" was used for points beyond the ideal boundary, which are now called hyper- 
ideal points. If we consider Delaunay triangulations up to similarity and hyperbolic 
polyhedra up to isometry, then the construction above establishes a 1-to-l cor- 
respondence between planar Delaunay triangulations and convex ideal polyhedra 
with one marked vertex. (Essentially the same construction, but represented in the 
projective model of hyperbolic space with a paraboloid as the absolute quadric, is 
known in Discrete Geometry as the "convex hull construction".) 

Weighted Delaunay triangulations are a well known generalization of Delaunay 
triangulations [15], where the sites are not points but circles (vertex-circles). We 
consider only the case where the closed disks bounded by the vertex-circles are pair- 
wise disjoint. Instead of an empty circumcircle, to each triangle there corresponds 
a circle (face-circle) which intersects the adjacent vertex-circles orthogonally and 
which intersects no vertex-circle more than orthogonally (see Figure 2). Weighted 
Delaunay triangulations with non-intersecting vertex circles correspond to hyper- 
bolic polyhedra with hyperideal vertices but with edges still intersecting hyperbolic 
space. Hyperideal points are not represented as such in the Poincarc half-space 
model. In the projective model, they are simply represented by the points outside 
the absolute quadric. The plane that is polar to such a point (with respect to 
the absolute quadric) intersects the absolute quadric, hence it represents a hyper- 
bolic plane. Thus, hyperideal points are in 1-to-l correspondence with hyperbolic 
planes. Two hyperbolic planes intersect orthogonally iff one is incident with the 
hyperideal point polar to the other. The correspondence between Delaunay trian- 
gulations and hyperbolic polyhedra extends to weighted Delaunay triangulations: 
A weighted Delaunay triangulation with non-intersecting vertex circles corresponds 
to a convex hyperbolic polyhedron with one marked ideal vertex and all other ver- 
tices hyperideal, and with edges intersecting hyperbolic space. The problem we will 
consider is to construct such polyhedra with prescribed combinatorial type and 
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dihedral angles. More generally, we will also consider weighted Delaunay triangula- 
tions in piecewise flat surfaces with cone singularities with prescribed cone angle at 
the vertices. These correspond to certain non-compact hyperbolic cone manifolds 
with polyhedral boundary where the lines of curvature connect one marked ideal 
vertex with the other hypcrideal vertices. 

1.3. Related work. For a comprehensive bibliography on circle packings and circle 
patterns we refer to Stephenson's monograph [34]. Here, we can only attempt to 
briefly discuss some of the the most important and most closely related results. 

Let C be a ccUulation of the 2-sphere and suppose there is a weight 9^ G (0, tt) 
attached to each edge e. Does there exist a hyperbolic polyhedron that is combina- 
torially equivalent to C and whose exterior dihedral angles are the weights 9e ; and if 
so, is it unique? A complete answer to this question is unknown. Andreev gave an 
answer for compact polyhedra with non-obtuse dihedral angles [2] (see also [29]), 
and he extended his result to non-obtuse angled polyhedra with finite volume, some 
or all vertices of which may be in the sphere at infinity [3] . An analogous existence 
and uniqueness theorem for circle patterns in surfaces of non-positive Euler char- 
acteristic is due to Thurston [35]. The intersection angles have to be non-obtuse 
but this theorem also allows the circle pattern equivalent of hypcrideal vertices. 
Chow & Luo [10] gave a proof which is inspired by work on the Ricci flow on sur- 
faces. They also show that there is a variational principle for this circle pattern 
problem, but only the derivatives of the functional are known. 

Rivin classified convex ideal polyhedra without any restriction to non-obtuse 
dihedral angles: 

Theorem 2 (Rivin [27]). There exists an ideal polyhedron that is comhinatorially 
equivalent to a cellulation C of and that has prescribed exterior dihedral angles 
Of, G (0, tt) iff for each cycle 7 in the 1-skeleton of the dual cellulation C* the 
inequality 



e* G7 

holds and equality holds iff 7 is the boundary of a face of C* . If it exists, the 
polyhedron is unique. 

Bao and Bonahon generalized Rivin's result for polyhedra with all vertices on 
or beyond the sphere at infinity: 

Theorem 3 (Bao & Bonahon [4]). There exists a polyhedron with ideal and hy- 
pcrideal vertices and with edges intersecting hyperbolic space that is comhinatori- 
ally equivalent to a cellulation C of and has prescribed exterior dihedral angles 
Oe G (0, tt) iff the following conditions hold: 

(i) For each cycle 7 in the 1-skeleton of the dual cellulation C* the inequality (1) 
holds and equality may hold only if 7 is the boundary of a face of C* . 

(ii) For each simple path 7 in the 1-skeleton of the dual cellulation C* that joins 
two different vertices of a face v* of C* and that is not contained in the boundary 



If it exists, the polyhedron is unique. 

Andreev, Rivin, and Bao & Bonahon obtain their results by employing variants 
the method of continuity (or deformation method) that was pioneered by Alexan- 
drov [1]. Schlenker gave a different proof of Bao & Bonahon's Theorem [31]. 

A variational approach to construct ideal hyperbolic polyhedra and, more gen- 
erally, Delaunay triangulations of piecewise flat surfaces was also provided by 
Rivin [26]. The basic idea is to build an ideal polyhedron by gluing together ideal 
tetrahedra, or equivalently, to build a Delaunay triangulation by gluing together 
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triangles. The angles of the triangles are considered as variables. They have to 
satisfy simple linear equality and inequality constraints: They have to be positive 
and the three angles in each triangle have to sum to tt. The angles around a vertex 
have to sum to 27r (more generally, some specified cone angle). Finally, to get the 
right circle intersection angles, the angles opposite an edge e have to sum to tt — Re- 
using Colin de Verdicre's terminology [12], we call an assignment of values to the 
angle variables that satisfies these constraints a coherent angle system. A coherent 
angle system does in general not represent a Delaunay triangulation because for 
the triangles to fit together, further non-linear conditions on the angles have to be 
satisfied. Nevertheless, the following theorem holds. 

Theorem 4 (Rivin [26]). Let S be a cellulation of a surface, let an intersection 
angle 9^ G [0, tt) be assigned to each edge e and a cone angle S.„ be assigned to each 
vertex v. There exists a Delaunay triangulation of a piecewise flat surface that is 
combinatorially equivalent to S and has intersection angles 9 and cone angles S if 
and only if the above constraints on the angle variables are feasible, i. e. if a coherent 
angle system exists. In that case, the Delaunay triangulation is unique up to scale. 

Note that the necessary and sufficient conditions of Theorem 2, involving inequal- 
ities of sums of intersection angles over paths, are very different from the conditions 
of Theorem 4, involving the existence of a coherent angle system. There is also a 
third type of conditions for the existence of a Delaunay triangulation of a piece- 
wise flat surface with prescribed intersection angles and cone angles that was flrst 
obtained by Bowditch [8] . It is by no means a triviality to directly derive one type 
of conditions from another type [28] [6]. A variant of Rivin's variational approach 
for Delaunay decompositions of hyperbolic surfaces was developed by Leibon [19]. 
In this article, we extend Rivin's variational approach to euclidean weighted De- 
launay triangulations with non-intersecting vertex circles. The main Theorem 1 is 
of the type "a weighted Delaunay triangulation exists uniquely iff a coherent angle 
system exists." The scope of Theorem 1 has non-empty intersection with Bao & 
Bonahon's Theorem 3. Both cover hyperbolic polyhedra with hyperideal vertices 
and precisely one ideal vertex. But the conditions are of a different type and not 
obviously equivalent. Theorem 1 also covers weighted Delaunay triangulations in 
flat tori and, more generally, in piecewise flat surfaces, possibly with boundary. The 
variational method is better suited for numerical computation. 

Recently, Schlenkcr has treated weighted Delaunay triangulations in piecewise 
fiat and piecewise hyperbolic surfaces using a deformation method [32] . He obtains 
an existence and uniqueness theorem [32, Theorem 1.4] with the same scope as 
Theorem 1, but the conditions for existence are in terms of angle sums over paths 
like in Theorem 3. This seems to be the first time that this type of conditions 
was obtained for circle patterns with cone singularities. It would be interesting to 
show directly that the conditions of his theorem are equivalent to the conditions of 
Theorem 1. 

Circle patterns have been applied to map 3D triangle meshes quasi-conformally 
to the plane. To improve these methods was an important motivation for this work. 
The group around Stephenson first used circle packings (circles touching without 
overlap) to construct planar maps of the human cerebellum [16]. This original 
method only takes the combinatorics of the input mesh into account. A later ver- 
sion uses so called inversive distance packings [9]. (The inversive distance of two 
non-intersecting circles is the cosh of the hyperbolic distance of the two planes they 
represent.) Inversive distance packings are similar to the weighted Delaunay trian- 
gulations considered here except that the inversive distances of the vertex-circles 
are prescribed instead of the intersection angles of the face-circles, and there is no 
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Delaunay criterion. Unfortunately, no existence and uniqueness theorem for inver- 
sive distance packings is known. Kharevych et al. [18] proceed along a different 
path. They first read off the intersection angles between circumcircles of the 3D 
triangle mesh. Then they construct a planar Delaunay triangulation with intersec- 
tion angles as close to the measured angles as possible. To construct the Delaunay 
triangulation they use a variational principle by Bobenko & Springborn [6], which 
is related to Rivin's via a Legendre transformation. It has the advantage that the 
variables arc (logarithmic) circle radii which arc not subject to any constraints. 
However, the prior angle-adjustment is achieved by solving a quadratic program- 
ming problem, which is the most complicated and computationally most expensive 
stage of the algorithm. One may hope that using weighted Delaunay triangulations 
will provide a way to escape the tight constraints that have to be satisfied by the 
intersection angles of a Delaunay triangulation without giving up all mathematical 
certainty regarding existence and uniqueness. Another interesting question is this: 
Can one formulate a dual variational principle for weighted Delaunay triangula- 
tions, with an explicit formula for the functional, where the variables are circle 
radii and inversive distances? 
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2. Euclidean hyperideal circle patterns 

2.1. Basic definitions. A surface is a two-dimensional manifold, possibly with 
boundary. A triangulated surface T (or triangulation for short) is a two-dimensio- 
nal CW complex whose total space is a surface S and which has the property that 
for each two-cell attaching map a : ^ S the set a~^{V) contains three points, 
where V is the vertex set (zero-skeleton) of the CW complex. 

This definition allows non-regular triangulations. A cell complex is called regular 
if the cell attaching homomorphisms embed the closed cells. A cell complex is 
called strongly regular if it is regular and the intersection of two closed cells is 
empty or a closed cell. The usual definition of simplicial complexes implies that 
they are strongly regular. Throughout this paper we assume all triangulations to 
be regular, but only to simplify notation. We will label vertices by i, j, fc, . . . 
and denote edges by pairs ij and triangles by triples ijk. However, this regularity 
assumption is not essential for the material presented here; everything holds for 
non-regular triangulations as well. 

We denote the set of vertices, edges, and triangles of a triangulation T by V, E, 
and T, respectively. Throughout this paper, all triangulations are assumed to be 
finite. 

A triangulated piecewise fiat surface [T , d) is a triangulated surface T equip- 
ped with a metric d such that for each two-cell attaching map a : —^ S the 
closed disk B^ equipped with the pulled back metric a*d is a euclidean triangle, 
and (T maps the vertices of this triangle to vertices of the CW complex. In other 
words, a triangulated piecewise flat surface is a surface obtained by glueing together 
euclidean triangles along their sides. (Of course, sides that are identified by glueing 
must have the same length.) The metric d is flat except in the vertices of the 



VARIATIONAL PRINCIPLE FOR WEIGHTED DELAUNAY TRIANGULATIONS 7 

triangulation where it may have cone-hke singularities. The cone angle at a vertex 
i is the sum of all triangle angles incident at i. If the cone angle at a vertex is 27r 
then the metric is flat there. A triangulated piecewise flat surface is determined by 
the triangulation T and the function I : E ^ Il>o that maps each edge ij £ E to 
its length kj. For each triangle ijk S T, the lengths kj, Ijk, hi satisfy the triangle 
inequalities. Conversely, a triangulation T and a function I : E ^ ]R>o that satisfies 
the triangle inequalities for each triangle determines a triangulated piecewise flat 
surface. 

A euclidean hyperideal circle pattern is a triangulated piecewise flat surface to- 
gether with a function r :V ^ lR>o with the following two properties. 

(i) For each edge ij £ E, r^ + rj < lij, where Ijj is the length of the edge. 

Let ijk £ T he a. triangle of the triangulation T. If we draw a triangle with sides 
hj, Ijk and Iki in the euclidean plane and circles with radii r^, rj and around 
the vertices, then the property (i) simply says that these circles do not touch or 
intersect. Consequently there exists a unique fourth circle that intersects all three 
circles orthogonally. The second condition concerns these orthogonally intersecting 
circles. 

(ii) Let ij G E he an interior edge. Let ijk and jil be the adjacent triangles on 
either side. (These may actually be one and the same triangle if the triangu- 
lation is not regular.) Draw two abutting triangles with the same side lengths 
in the euclidean plane, and draw circles with radii r^, rj, and n around the 
vertices. Then the orthogonal circle through the vertex-circles of one triangle 
intersects the fourth vertex-circle either not at all or at an angle that is less 
than ^. 

In other words, a euclidean hyperideal circle pattern is a weighted Delaunay 
triangulation with non-intersecting vertex-circles in a piecewise flat surface. (Note 
that condition (ii) invokes an edge-local Delaunay condition. This raises the ques- 
tion whether it is also true for piecewise flat surfaces that the local condition implies 
the global condition that no face-circle intersects any vertex-circle more than or- 
thogonally. Such questions shall not be treated here. We refer to Bobenko & Spring- 
born [7] and Bobenko & Izmestiev [5] for a more thorough treatment of Delaunay 
triangulations and weighted Delaunay triangulations in piecewise flat surfaces.) 

Just as a triangulated piecewise flat surface is obtained by glueing together 
euclidean triangles along the edges, a euclidean hyperbolic circle pattern is obtained 
by putting together triangles which are decorated by circles as shown in Figure 3. 

2.2. Interpretation as hyperbolic polyhedra. By the construction explained 
in Section 1.2, the decorated triangle shown in Figure 3 may be interpreted as a 
tetrahedron in three-dimensional hyperbolic space with one vertex on the sphere at 
infinity and three vertices beyond that sphere, as shown in Figure 4. The sides of 
the triangle and the face-circle that intersects the three vertex-circles orthogonally 
correspond to hyperbolic planes that bound a tetrahedron with one ideal and three 
hyperideal vertices, which are represented by the vertex-circles. 

If we can put together the decorated triangles to form a circle pattern, we can also 
glue the corresponding hyperbolic tctrahcdra to form a hyperbolic cone manifold 
with polyhedral boundary. Such a glueing will identify the ideal vertices of all 
tetrahcdra. The resulting point will either be an ideal vertex if T has boundary 
or a cusp of the hyperbolic manifold if T is closed. The hyperideal vertices of 
the tetrahedra will be identified in groups corresponding to vertices of T to form 
hyperideal vertices of the polyhedral boundary. There will in general be cone lines 
running from the cusp or ideal vertex to the hyperideal vertices. If the circle pattern 
has no curvature at the vertices, we obtain a non-compact hyperbolic manifold with 
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Figure 3. Local geometry at a triangle. 
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Figure 4. A tetrahedron in hyperbolic space (shown in the projective model) with 
one ideal vertex and three hyperideal vertices. The tetrahedron is truncated by the 
polar planes of the hyperideal vertices. 

polyhedral boundary. If in addition the triangulation T is topologically a disk, we 
obtain a convex hyperbolic polyhedron with one ideal vertex all others hyperideal. 

2.3. The circle pattern problem. From a euclidean hyperideal circle pattern 
one can read off the following data: 

• A triangulated surface T. 

• For each vertex i g y, the sum S (0, cxd) of incident triangle angles. For an 
interior vertex, this is the cone angle at z, and Si = 27r if the circle pattern is flat 
at i. For a boundary vertex, is the interior angle of the polygonal boundary 
at i. 

• For each interior edge ij £ i?, the intersection angle 6*^ G [0, tt) of the orthogonal 
circles corresponding to the two adjacent triangles zjfc and jil. An intersection 
angle 6*^ = means that the orthogonal circles of triangles ijfc and ji^ coincide. 
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• For each boundary edge ij G the intersection angle, also denoted by 9ij, of the 
orthogonal circle corresponding to the adjacent triangle ijk with the line segment 
containing the edge ij. The range of these intersection angles at boundary edges 
is 9ij e (0,7r). 

We consider the following circle pattern problem: Given an abstract triangulation 
T and angle data E : V (0, oo), 9 : E ^ [0,Tr) find a corresponding euclidean 
hyperideal circle pattern. 

2.4. Local geometry at a triangle. Consider a geometric figure consisting of a 
euclidean triangle with three non-touching and non-intersecting circles centered at 
the vertices and a fourth circle intersecting the other three orthogonally. Let ai2, 
0^23, Q;3i, 7i, 72, 73 be the angles shown in Figure 3. They are positive, 

aij > 0, 7j > 0, (2) 

satisfy the angle sum equation 

71 + 72 + 73 = 71", (3) 

and the inequalities 

71 + ai2 + 0:31 < 

72 + ^23 + "12 < Ti", (4) 

73 + 0:31 + Ck23 < TT. 

Let 

A = {(ai2,a23,a3i,7i7 72,73) £ satisfying (2), (3), and (4)} . 

Conversely, if (q;i2, a23, "31, 7i, 72, 73) & A, then there exists one such figure 
with these angles, and only one up to similarity. Indeed, the construction of such a 
figure is simple: Draw any circle in the plane. (This fixes the scale and translational 
degrees of freedom.) Then draw a line intersecting it at the angle ai2. (This fixes the 
remaining rotational degree of freedom.) Then draw the other two lines intersecting 
the circle and the first line at the prescribed angles. The inequalities (4) ensure 
that the lines intersect outside the face-circle and that the orthogonal vertex-circles 
do not intersect each other. 

2.5. Coherent angle systems. Let T be a triangulation with triangle set T, edge 
set E, and vertex set V. We label the coordinates of points in R^l-^l by a*^ , 7* where 
t €T is a triangle and i,jQt are vertices of t. Wc fix this labeling once and for all. 
Let S : 1/ — > (0, 00) be a function on the vertices and 6' : i? [0, tt) be a function 
on the edges. The space of coherent angle systems -4(T, 9, S) is the set of all points 
(a,7) e E.^l'^l such that 

• for each triangle t ~ ijk G T 

ft t t t t t\ ^ \ 
(ay,ajfc,afe,,7i,7j,7fc) £ A, 

• for each interior edge ij € E 

a*^ =7r-%, (5) 

where t,t' are the adjacent triangles on either side of edge ij, 

• for each boundary edge ij G E, 9ij > and 

a^=7r-%, (6) 

where t is the triangle incident with edge ij, 

• for each vertex i £ V 

E 7* = S.. (7) 

tGT:iet 
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A coherent angle system is an element of A{T, 6, S). If A{T, 6, ^) is not empty, 
then the closure A{T,9,'E.) is a compact polytope in R^l^l, and A{T,9,'E.) is its 
relative interior. 

2.6. Existence and uniqueness. The main Theorem 1 reduces the question of 
existence and uniqueness of a solution of the circle pattern problem of Section 2.3 to 
a linear feasibility problem. The "only if" part of the theorem — if a circle pattern 
exists then a coherent angle system exists — is trivial, since one can simply read off a 
coherent angle system from a euclidean hyperideal circle pattern. The "if" part — if 
a coherent angle system exists then a circle pattern exists — is not trivial, and the 
rest of this paper is devoted to proving it. It is not true that for each coherent angle 
system there exists a euclidean hyperideal circle pattern with these angles. While it 
is true that a coherent angle system determines up to similarity a geometric figure 
as shown in Figure 3 for each triangle of the triangulation, these figures cannot 
in general be put together to form a circle pattern. The relative scale of the two 
figures corresponding to neighboring triangles is determined by the condition that 
the triangle edges to be glued together must have the same length. But the two 
pairs of corresponding vertex circles in the two figures will in general not have 
matching radii. Even if we disregarded the vertex circles, it is in general not be 
possible to choose consistently a scale for each triangle such that the corresponding 
sides of triangles that are to be glued together have the same length. 

3. A VARIATIONAL PRINCIPLE 

The space ^(T, 0,5) C R^'^l of coherent angle systems is defined by linear equa- 
tions and inequalities. For a coherent angle system to describe a solution for the 
circle pattern problem, it has to satisfy in addition certain non-linear equations, 
which guarantee that the triangle figures can be put together to form a circle pat- 
tern. It turns out that these non-linear compatibility conditions are the conditions 
for a critical point of the function F : R^'-^l R defined below under variations 
in yl(T, 0, S) C R^l'^l. This is the content of Lemma 1. Together with Lemmas 2 
and 3 this provides a proof of Theorem 1. 

3.1. The functional. Label R^'-^l as in Section 2.5 and define 

Ft : R*^l^l ^ R 



by 



where 



Ft^ ^K-'"5fc'"L,7-,7*,7D, (8) 

t=ijk£T 
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-a23— Q^3i- 


-73 



(9) 



JI(a;) = - /" log|2sine|d^ 
Jo 



is Milnor's Lobachevsky function [22], [24]. (This is up to scale the same as Clausen's 
integral Cl2(a;) = 2JI(|); see Clausen [11], Lewin [20].) The function JI is tt- 
periodic, continuous, and odd. It is smooth everywhere except at integer multiples 
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Figure 5. Milnor's Lobachevsky function, y = JI(a;). 



of TT where its graph has a vertical tangent; see Figure 5. We will simply write F 
for Fq- when the triangulation T can be inferred from the context. 

3.2. The main lemmas. The following three lemmas imply Theorem 1. By 
Lemma 1, the critical points of F in ^(T, 6',S) correspond to the solutions of 
the circle pattern problem. Lemma 2 implies the uniqueness claim of Theorem 1. 
The existence claim follows from Lemma 3. 

Lemma 1. A coherent angle system p G ■4(T, 0, S) is a critical point of F under 
variations in A{T, 9, S) if and only if the decorated triangles fit together. 

Lemma 2. The function F is strictly concave on A(T, 6, S). 

Lemma 3. //yl(T, 0,S) is non-empty, then the restriction of F to the closure 
yl(T, 0,S) attains its maximum in ^(T, 0,S). 



The key ingredient to this proof is Schlafli's formula for the derivative of the volume 
of a hyperbolic polyhedron when it is deformed in such a way that its combinatorial 
type is preserved; see Milnor [25], [30]. 

Theorem (Schlafli's differential volume formula). The differential of the volume 
function V on the space of 3-dimensional hyperbolic polyhedra of a fixed combina- 
torial type is 



ij 

where the sum is taken over the edges ij, and Oij, Ui are the length and interior 
dihedral angle at edge ij. 

Hyperbolic polyhedra with some or all vertices on the sphere at infinity still have 
finite volume. Milnor notes (in the concluding remarks of [25]) that Equation (10) 
remains true for such polyhedra, under the following modification which is necessary 
because the edges incident with an ideal vertex are of course infinitely long. Choose 
arbitrary horospheres centered at the ideal vertices, and for edges ij incident with 
an ideal vertex let be the length of the edge truncated at the horosphere(s) 
centered the ideal endpoint(s). One should add that only deformations that leave 
the ideal vertices on the infinite sphere are considered. The angle sum at an ideal 
vertex remains constant under such a deformation, so ^ daij ~ 0. If one chooses a 
different horosphere at an ideal vertex i, the resulting truncated edge lengths of the 
incident edges differ by the same additive constant: becomes +Ci. Hence the 
right hand side of Equation (10) does not depend on the choice of horospheres. A 



4. Proof of Lemma 1 




(10) 
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proof for this extension of Schlafli's differential volume formula to polyhedra with 
ideal vertices is contained in [33] (Lemma 4.1). 

Next we extend Schlafli's differential volume formula to polyhedra that have 
vertices beyond the sphere at infinity, but all of whose edges still intersect hyperbolic 
space. If we truncate the hypcrideal vertices at the dual hyperbolic planes, we 
obtain a polyhedron with finite volume. (It may still have ideal vertices.) 

Definition. The truncated volume of a hyperbolic polyhedron with vertices beyond 
infinity is the finite volume of the corresponding truncated polyhedron, truncated 
at the hyperbolic planes dual to the hypcrideal vertices. 

Equation (10) remains true for polyhedra with hyperideal vertices if we let V 
be the truncated volume and a-y be the edge lengths truncated at the dual planes 
of hyperideal vertices and at horospheres centered at ideal vertices. Indeed, if we 
apply Schlaflis differential volume formula to the truncated polyhedron, the edges 
introduced by the truncation (those between original faces and truncation planes) 
have dihedral angle 7r/2, and this angle is constant during any deformation. So for 
these edges daij = 0. Thus, only terms involving the original edges remain in (10). 

Lemma 4. The truncated volume of the tetrahedron with one ideal and three hyper- 
ideal vertices with angles as shown in Figure 4 is equal to V{ai2, 023, asi, 7i, 72, 73) 
as defined by Equation (9). 

We prove this lemma in Section 7. 

Remark. Exactly the same formula holds for tetrahedra with one ideal and three 
finite vertices [36] , [37] . The only difference is that the dihedral angles for such poly- 
hedra satisfy the opposite inequalities instead of (4). This seems to be a common 
phenomenon: that the same volume formula holds regardless of whether certain 
vertices are finite or hyperideal. In the case of orthoschemes (the generalization of 
triply orthogonal tetrahedra to arbitrary dimension) this was observed by Keller- 
hals [17]. Section 7 contains some more examples. It would obviously be useful to 
have this observation cast into a theorem. The author is not aware that this has 
been done. 

From Schlafli's differential volume formula and Lemma 4 we obtain the following. 

Lemma 5. // we choose a horosphere centered at the ideal vertex of the truncated 
hyperbolic tetrahedron shown in Figure 4 then 



where Oij is the length of the edge between the hyperideal vertices i and j truncated 
at the polar planes, and 



where Ui is the length of the edge from the ideal vertex to the hyperideal vertex i, 
truncated at the horosphere and the dual plane. 

Equations (11) and (12) provide formulas for the edge lengths of the truncated 
hyperbolic polyhedron in terms of the dihedral angles and 7^. Because the 
choice of the truncating horosphere at the ideal vertex is arbitrary, the lengths Oi 
are only determined up to an additive constant. 

The hyperbolic lengths and ai arc related to the cuclidcan lengths Uj and 
the radii r.i (see Figure 3). The radii are proportional to e^"*, i.e. 




(12) 




(13) 
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and 

=rl+ri + 2rir2 coshaia. (14) 
(These relations arc obtained by straightforward calculation in the Poincarc half- 
space model. The quantity coshai2 is the inversive distance of two circles in the 
plane.) Together, Equations (11)-(14) provide formulas for the radii and the 
euclidean edge lengths kj in terms of the angles and -ft . They determine the 
and lij up to a common factor, in agreement with the fact that the angles determine 
the decorated triangle of Figure 3 up to similarity. 

Now let {a, 7) S A{T, 9, 5) be a coherent angle system. For each triangle t G T 
the angles (a-j , a*^., a^j, 7*, 7*, 7^) determine a decorated triangle up to similarity. 
These fit together to form a hyperideal circle pattern iff they can be scaled consis- 
tently; this means iff a radius can be assigned to each vertex i and a length lij to 
each edge such that the relations (11)-(14) hold for each triangle. Equivalently, the 
corresponding hyperbolic tetrahedra fit together iff the horospheres at the infinite 
vertices can be chosen consistently; this means iff an G E, can be assigned to 
each vertex i and an S E, to each edge ij such that the following holds: 

1ft £T and i,j G t, then 

„ dV 



'da*- - " 



(15) 



and 



3 

where the derivatives are evaluated at (a*j , ctj^, ' 7i ' 7j' 7a:)- Now since 
— p {a,^ , a^-fc , a^, , 7, , 7,: , 7^) = — ^ {a, 7) 

ij ij 

and 

^(ay,a^fc,afc,;,7,;,7j,7fc) = ^("'^)' 

(see Equation (8)), Equations (15) and (16) are equivalent to 

OF 

- 2^r;(-(a:7) = (17) 

ij 

and 

fdF OF \ 

"Ha^^"'^^" ^^^^ 

Clearly, the Equations (17) for the a.y are compatible iff the following condition 
holds: 



(i) For each interior edge ij, 

/ d d 



r)F(a,7) = 0, (19) 



\da\j da]^ 

where t and t' are the triangles adjacent with ij . 

Equations (18) for the ai are compatible iff the condition holds: 

(ii) If iotiiit2i3 ■ ■ ■ tnin is any finite sequence of alternatingly vertices and trian- 
gles that starts and ends with the same vertex io = i„, and that has the property 
that each contains preceding vertex im-i and the following vertex i„i, then 

These conditions (i) and (ii) are the non-linear compatibility conditions that a 
coherent has to satisfy to represent a solution to the circle pattern problem. It 
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remains to show that they are also the conditions for a critical point of F under 
variations in .4(T, 0, S). This is achieved by the following lemma, which concludes 
the proof of Lemma 1 . 

Lemma 6. The tangent space to A{T , 6, S) is spanned by the tangent vectors 

d d 



(one for each interior edge) and 



y.{A — (22) 

(one for each cycle 10^1*1^2^3 • ■ - tnin) that appear in conditions (i) and (ii) above. 

Proof. The space A{T , 6, S) is defined by strict inequalities and Equations (3), (5), 
(6), and (7). The equations for a tangent vector are therefore: 

(a) For each boundary edge ij: da\j = 0. 

(b) For each interior edge ij ^ tDt': dajj + da*^^ = 0. 

(c) For each vertex i: c?7* = 0. 

(d) For each triangle t = ijk: djf + ^7* + djl ~ 0. 

Since each equation involves either alphas or gammas but not both, the tangent 
space is the direct sum of an a-subspace and a 7-subspace. The a-subspace is 
clearly spanned by the tangent vectors (21). To see that the 7-subspace is spanned 
by the vectors (22), let G be the graph with vertex-set V ^ V U T and edge-set 

£ = {{i,t} eV\i eV,t eT,i et}. 

The edges {i,t} of G are in one-to-one correspondence with the tangent vectors 
. This gives rise to a linear isomorphism between the space of edge-chains of G 

(over R) and the space spanned by the vectors . Equations (c) and (d) are then 
simply the equations for the cycle-space of G. □ 

5. Proof of Lemma 2 

We are going to show that each of the terms y(a-, a*, a^, 7*, 7*, 7^) in Equation (8) 
is concave. To this end, we split the 15-term sum in Equation (9) which defines 
V into five parts; see Equation (25). Each part represents the volume of an ideal 
tetrahedron (Theorem 5), which is known to be concave (Lemma 7). 

5.1. The volume of an ideal tetrahedron. Consider a hyperbolic tetrahedron 
with all four vertices on the sphere at infinity. For each vertex the sum of the 
interior dihedral angles at the adjacent edges is tt. This implies that the angles at 
opposite edges are equal [22] [24] . An ideal tetrahedron is therefore determined by 
three angles in the set 

Ao = {(71,72,73) e K.^ I 7» > 0, 71 +72 + 73 = tt}. (23) 

Theorem 5 (Milnor [22] [24]). The hyperbolic volume of an ideal tetrahedron with 
dihedral angles (71,72,73) G Ao is 

Vb(7i,72,73) = JI(7i) + JI(72) + JI(73)- (24) 

Lemma 7 (Rivin [26]). The volume function Vq is strictly concave on Aq. 

The proof is straight forward. For the reader's convenience, we repeat it here. 
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Proof of Lemma 7. Let 

/(a, [3) = Vo{a,f3,Tr-a-P) = Jl{a) + JI(/3) - Jl{a + (3) 

and assume that a > 0, /3 > 0, a + /3 < tt. Since J\" {x) = — cotx, the Hessian 
matrix of / is 

- cot a + cot(a + /?) col{a + (3) 



Hess / 



cot(a + /3) — cot (3 + cot{a + j3) 



A short calculation shows that the determinant of Hess(/) is 1. The matrix is 
therefore either positive definite or negative definite. But since the cotangent is a 
strictly decreasing function on (0,7r), the diagonal elements are negative. Hence 
Hcss(/) is negative definite and /(a,/?) is strictly concave. □ 

5.2. Five ideal tetrahedra. Equation (9) for the truncated volume V{ai2, a23, Q^ai, 7i, 72, 73) 
can be rewritten as 



2ViT) = ^ ( JI(7,:) + Jl(7.:) + -n(7-') + + 



where 



i=l 



7^ = 2 ^\ 72 = 2 ' 73 

// _ 7r-a3i+ai2— 71 _ 7r— 012+023-72 

7l — 2 ' 12 — 2 ' '3 



2 




7r — 031 +012 


-71 


2 




''r + Ct31+C<12 


-71 


2 




IT— Q31 — ai2 


-71 


2 



7r+Qi2 — ^23 


-72 


2 




TT — ai2+a23 


-72 


2 




7r + Ql2+a23 


-72 


2 




7r — ai2 — CK23 


-72 



7r- 


-f 023—^31 


-73 




2 




7r 


— 023 + ^31 


-73 




2 




7r- 


-f a23+a3l 


-73 




2 




7r- 


— Q23— "31 


-73 



Ml — ^ — 2 ' — 2"^ , M3 

11 — 031—012—71 TT- Q!12— 023— 72 ,. 

^1 2 ' 2 2 ' 3 

The following observation is both very simple and crucial for this proof (and also 
for the proof of Lemma 3 in Section 6): If 

(ai2, 023,031, 71,72, 73) e A, 

then 

(71,72,73) gAo, (7;',72,73)eAo 

and also 

(71,^1,1^1) e Ao, (72,^2,1^2) e Ao, (73,^3,1^3) e Aq. 
Thus, 2V is the sum of the volumes of five ideal tetrahedra: 

21/(012,023,031,71,72,73) = K)(7i,72,73) + ^0(7", 72, 73) 

+ V'o(7i,Mi,'^i) + Vb(72,M2,i^2) + Vb(73,/i3, i^3)- (25) 

Since each of the five terms is concave by Lemma 7, 2V is concave and so is F. 
This completes the proof of Lemma 2. 

Remark. We have no geometric explanation why 2V is the sum of the volumes of 
five ideal tetrahedra, although this may well be a consequence of Doyle & Lei- 
bon's "23040 symmetries of hyperbolic tetrahedra" [14]. Equation (25) is not 
the only way to write 2V as a sum of five tetrahedra: For example, because 
(ai,Q;2,a3,7i,72,73) S A also implies 

(71,72,73) e Ao, (7i,/i2,i^3) e Ao, (72,^3,1^1) e Aq, (73,/ii,:/2) e Aq, 

one has as well 

2F(ai2,023,03i,7i,72,73) = Vb(7i,72,73) + Vo(7i', 72 , 73 ) 

+ ^o(7l,M2,J^3) + Vo(72,/"3,l^l) + Vo(73,/ii, J/2). 
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6. Proof of Lemma 3 



To prove Lemma 3 we have to show the following: 

Claim. Suppose A{T, 6*, S) ^ and let p G A{T,e,E) \ A{T, 9, S). Then there is 
aqe A{T, 9, S) with F{q) > F{p). 

In the following Section, we will analyze the behavior of the volume function V 
as the dihedral angles approach the relative boundary of the domain. In Section 6.2 
we will use this analysis to prove the Claim. 

6.1. Behavior of the volume function at the boundary of the domain. We 

will first recollect Rivin's analysis [26] of the behavior of the volume Vo(p) of an 
ideal tetrahedron (see Section 5.1) as p approaches the relative boundary Aq \ Aq of 
the domain Aq. Then we will use this for the corresponding analysis of the volume 
function y of a tetrahedron with one ideal an three hyperideal vertices. Here we will 
again make essential use of the decomposition into five ideal tetrahedra described 
in Section 5.2. 

Boundary behavior of Vq . First consider the set Aq of dihedral angles of ideal tetra- 
hedra (see Equation (23) in Section 5.1). Its closure Aq is the 2-simplex in E,^ that 
is spanned by the points (7r,0, 0), (0, 7r,0), (0,0, tt). The points in the relative 
boundary Aq \ Aq correspond to ideal tetrahedra that have degenerated to planar 
figures. 

Definition. We call a point (71,72,73) G Ao \ Ao mildly degenerate iff (71,72,73) 
is some permutation of (0, /3, tt — /3) with < /3 < tt, i.e. iff p is contained in an 
open side of the boundary triangle. We call it hadly degenerate iff (71,72,73) is 
some permutation of (0, 0, tt), i.e. iff it is a vertex of the boundary triangle. 

It is easy to see that Vb vanishes on Aq \ Aq. But the speed with which V{q) 
tends to as g G Aq approaches the boundary is different depending on whether q 
approaches a mildly or a badly degenerate point: 

Lemma 8. Let p G Aq \ Ao, q G Aq. If p is mildly degenerate, then 



has a finite positive limit for t \ 0. 

Proof. The claim for mildly degenerate p follows directly from the fact that JI (x) is 
smooth except at integer multiples of tt, where the derivative J\\x) = — log |2 sinxj 
tends to +00. To prove the claim for badly degenerate p, let us assume without 
loss of generality that p = (0, 0, tt). Then 



Mm —Vo{{l - t)p + tq) ==+cx). 



(26) 



If p is badly degenerate then the t-derivative 



±Vo{{l ^t)p + tq) 



(27) 



(1 - t)p + tq = {ta, tb, TT -t{a + b)) 

for some a, 6 > 0, and 



Vo{{l - t)p + tq) = Jl{ta) + Jl{tb) - JI {t{a + b)). 



Hence for the t-derivative we obtain 




a+b 



> 0. 



□ 
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Boundary behavior of V . Now consider the set A of dihedral angles of tetrahedra 
with one ideal and three hyperideal vertices (see Section 2.4). Its closure A is the set 
of points (q!i2, a23, asi, 71, 72, 73) G K.^ that satisfy Equation (3) and the non-strict 
versions of the inequalities (2) and (4). A point is in the relative boundary A\ A if it 
satisfies with equality at least one of these the non-strict inequalities. In Section 5.2 
we described an affine map A (Aq)^ associating five ideal tetrahedra with each 
point in A. This extends to a map A — ^ (Aq)^. We classify the degenerate points 
p G A \ A according to whether any of the five corresponding ideal tetrahedra 
degenerate and the way in which they do: 

Definition. Let p G A \ A. We say that p is mildly degenerate iff at least one of 
the five corresponding ideal tetrahedra is mildly degenerate. We say that p is badly 
degenerate iff at least one of the five ideal tetrahedra is badly degenerate but none 
are mildly degenerate. We say that p is a-degenerate iff all five corresponding ideal 
tetrahedra are non-degenerate. 

Clearly, every p G A \ A is either mildly degenerate, badly degenerate or a- 
degenerate. By Lemma 10, there is only one badly degenerate p up to a permutation 
of the indices. The reason for calling the last type "a-degenerate" will be made 
clear by Lemma 11. The next lemma follows immediately from Lemma 8. 

Lemma 9. Let p e A \ A, g G A. If p is mildly degenerate, then 

\im^V{{l -t)p + tq) ^ +00. (28) 
If p is badly degenerate or a-degenerate then the t-derivative 

^V{{l-t)p + tq) (29) 

has a finite limit for t \, 0. 

Lemma 10. Suppose (ai2, 0123, 0131, 71, 72, 73) G A\ A is badly degenerate. Then 

li = oijk = TT, 7j = 7fc = aij = = 

for some permutation {i,j,k) of (1,2,3). (This implies that all five corresponding 
ideal tetrahedra are badly degenerate.) 

Proof. We consider separately all essentially different cases of one of the five ideal 
tetrahedra badly degenerating in some way, where we consider cases as essentially 
different if they do not only differ by a permutation of the indices. First, note that 
by summing (the non-strict versions of) the inequalities (4) and using Equation (3), 
one obtains 

ai2 -I- Q!23 + ^31 < 

Case 1: (71, 72-73) = (0,0, tt). Since 

I" = 273 - = Q!23 - ^31 - 73, 

we have 023 = tt and hence ai2 = a^i = 72 = 73 = and 71 tt. 

Case 2: (7i,/ii,j'i) = (7r,0, 0). First, 71 = tt implies 72 = 73 = 0^12 = q;3i = 0. 
Then 

I I I i\ (a ^ ~ "23 7r + a23 

(71,72,73) = (0, — ^ — , — ^ — ), 

and since by assumption this is not mildly degenerate, we have 023 = tt. 

Case 3: (71, /zi, vi) = (0, tt, 0). First, 71 = implies 72 + 73 + tt; /^i = tt implies 
c^3i +tti2 = TT and hence a23 = 0. Now the non-strict versions of the inequalities (4) 
imply Q!i2 -l- 72 = and a^i + 73 = tt. Hence 1^2 = ^^3 = 0. Since by assumption, 
neither (72,^2,^^2) nor (73, /is, 1^3) are mildly degenerate, either 72 = 0:31 = tt and 
73 = ai2 = or 72 = a^i = and 73 = ai2 = tt. 

(Non-)Case 4-' (71,^1,^^1) = (0,0, tt). This cannot happen because < ^. □ 
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Lemma 11. Suppose p = (ai2, 023, 031, 71, 72, 73) G A \ A is a-degenerate. Then 
all "fi > and the strict inequalities (4) are satisfied, but some Uij vanish. 

Proof. First, 7^ > and 7^ + a.y + aki < tt because 7^ = or 7.^ + aij + aki — tt 
would imply that (7^, ^i, Vi) is degenerate. But then some must vanish because 
otherwise p G A. □ 



6.2. Proof of the Claim. Suppose ^(T, 6*, S) ^0andletp= (a, 7) £A{T,6,'~)\ 
^(T, 6',S). We distinguish several cases. 

Case 1. There is at least one triangle t T such that (a*, 7*) is mildly degenerate. 
Let p = {a, 7) G ^(T, 0, S) be any coherent angle system. Then Lemma 9 implies 

lim - t)p + tp) = +00. 

*\o dt ^ ' 

It follows (by the mean value theorem) that if £ > is small enough, then > 
F{p) for g = (1 — s)p + ep. This completes the proof of the Claim under the 
assumption of Case 1. 

Case 2. There are no t G T with (a*, 7*) mildly degenerate or a-degenerate. 
This means all degenerate (a*, 7*) are badly degenerate. We will construct a 
p e A{T, 9, S) \A{T, 9, 5) which satisfies F{p) = F{p) and the conditions of Case 1. 
Since the Claim was proven for Case 1, it holds in Case 2 also. 

Suppose ti G T with (a*, 7*) badly degenerate. Let ij be the edge of ti with 
alj = TT. Equations (5) and (6) imply that 9ij = and therefore that edge ij 
is not a boundary edge. Let t2 & T he the triangle neighboring ti across edge 
ij. Again by Equation (5), a*^^ = 0, so t2 is also badly degenerate. By repeating 
this argument we construct a sequence ti,t2,t3, . . . of badly degenerate triangles, 
each one adjacent to the next. Since there are only finitely many triangles, this 
sequence must eventually loop back on itself. Let us reindex the triangles such 
that ti, t2, • ■ • , in = ii is such a loop of badly degenerate triangles. Define an angle 
system p = (a, 7) as follows (see Figure 6). If t G T is not contained in the loop 
of triangles, then let (a*, 7*) = (a*, 7*). If < = i,„ is contained in the loop, let 
ij = im n tm+i and jk tm-i n tm, hence 

a*j = 7fc 7^, "jfc = = Ji = Ij = 0. 

Let 







= TT — e. 




= It 






= 1] 


= 0, 



with e G (0, tt) arbitrary but the same for all triangles in the loop. We claim that 
p G A{T, 9, S) \ A{T, 9, S). Indeed, the triangles in the loop are still in A \ A, but 
now they are mildly degenerate instead of badly degenerate. Further, the sum of 
a-angles at edges has obviously not changed, so Equation (5) still holds. Finally, 
to see that the sum of 7-angles around each vertex has not changed, note that at 
a vertex i this angle sum is increased by e for each time that the loop of triangles 
enters the star of i and is decreased by e each time it leaves the star. So Equation (7) 
still holds. Also F{p) = F{p), because the truncated volumes for the triangles in 
the loop are before and after the deformation. Hence we have reduced Case 2 to 
Case 1. 



Case 3. There are no t G T with (a*, 7*) mildly degenerate, but there is at least one 
t gT such that (a*, 7*) is a-degenerate. Suppose t = ijk G T is a-degenerate and 
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Figure 6. How to modify the angles along a loop of badly degenerate triangles. We 
write TT— and 0+ as shorthand for "tt is replaced by tt — e" and for "0 is replaced by 
e" , respectively. The changes in a-angles are indicated along the edges, the changes 
in 7-angles are indicated in the corners of each triangle. 



alj = 0. Let t' = jil be the triangle on the other side of edge ij. Note that a*^ 
cannot vanish, because this would imply 9ij = tt — a*^ — a*j = tt. We distinguish 
two sub-cases. 

Case 3(a). t' is not badly degenerate. Then < a*^ < tt, and we can change the 
angle system p to the angle system p by setting a* ■ = e and a 



£ for some 



small e > 0, and keeping all other angles the same. We claim that F{p) > F{p) if 
e is small enough. Indeed, 

^ -y(a*,7*)=0, 



because V(a*, 7*) 



is even in the a- variables (see Equation (9)), and 



da 



2«i» < 0> 



where a*j is the length of the corresponding truncated edge (Lemma 5), which 
is positive because aji = only if aji = (see Equation (14) and Figure 3). So 
provided that e > is small enough, V {a* , + V {a^ ,7* ) increases if q;*^ increases 
by £ and a* ^ decreases by e. We have thus constructed a p € A{T , 9, S) without 
mildly degenerated triangles such that F(j)) > F{p) and the number of vanishing 
a-angles has decreased by one. The Claim follows by induction on the number of 
vanishing a-angles. 

Case 3(h). t' is badly degenerate We will construct a q € A{T,6,E) such that 
F{q) > F{p), but the construction is less straightforward than in the other cases. 



Note that in this case, 



= TT and hence On ~ 0. This means that the solution of 
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Figure 7. Left: If neighboring triangles share the same orthogonaUy intersecting 
circle, we can perform an edge flip. Right: How q{s,t) differs from q. 





Figure 8. The a-degenerate tetrahedron t 
triangles ti = {ilk), t2 = (jkl). 



(ijk) with a*- = is split into two 



the circle pattern problem (if it exists) is such that the same face-circle corresponds 
to both triangles t and t'; see Figure 7 (left). Equivalently, the tetrahedra corre- 
sponding to these triangles fit together to form a pyramid over a quadrilateral base. 
Thus, one could pose an equivalent circle pattern problem using the triangulation 
T that is obtained by flipping the edge ij. This is the basic motivating idea behind 
the construction we will now describe. 

First, we split the a-degenerate tetrahedron (a*, 7*) into two tetrahedra (a*^ , 7*"^) 
and (a*2,7*2) as shown in Figure 8. This introduces new angles 7^,^ 



satisfy 



and also 



which are uniquely determined by the old angles (a , 7 ). 



7fc% ri% 

They clearly 



ri 



''Ik 



"kl 



,*2 
^kl 



ll + ^Ik = ll 

The remaining angles in (q;*^,7*^) and (q;*^,7*^) are equal to the corresponding 
angles in (a*, 7*): 



7, 



a 



ki 



a 



ki ' 



a 



a 



0. 



The volumes satisfy 



y(a*,7*) = y(a*S7*0 + V'(a',7") 



(30) 
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The tetrahedra [a*^ , 7*^ ) and (a*^ , 7*^) are mildly degenerate, because they are not 
badly degenerate but 

1 I ^1 ^2 I ^2 1 ^2 

7; + a/fc + ot,] - TT, 7; = + + a^i = tt. 

Now let T be the triangulation obtained from T by flipping the edge ij, thus 
replacing it with an edge kl and replacing the triangles t, t' with triangles ii, 
(This edge flip can be performed even if the triangulation T is not regular. The only 
obstruction for an edge to be flippable is that it is adjacent to the same triangle on 
either side. But this is not the case here, because (a*, 7*) is a-degenerate, whereas 
(a* ,7' ) is badly degenerate.) Let p be the angle system with (0*^7*^), (a*^,7*^) 
as described above and all other angles the same as in p. Then 

pe^(r,0~,2), 

where Omn = Qmn for all edges ran of T except for kl and 9ki — Oij — 0. Because of 
Equation (30) and because the volume of the badly degenerate tetrahedron (a* , 7* ) 
vanishes, we have Fq-{p) = Ff{p). Now ^(T, 9,E) ^ (i) because from every coherent 
angle system in A{T, 0, S) one can easily construct a coherent angle system in 
A{'T, 9, S). Hence the reasoning of Case 1 above applies and there exists a 

g= (5,7) e^(f,^,S) with F^iq) > F^iP) ^ Frip). 

It remains to construct a. q G A{T,9,E) with Friq) > F^{q). If the tetrahedra 
(5;*% 7*1) and (5*^,7*2) fit together, this could be achieved by performing another 
edge flip as in Figure 7 (left). However, they will in general not fit together. We will 
therefore deform q to obtain a. q € A{T, 9, S) such that the triangles fit together 
and Ff{q) > F^{q). To this end, let q{s,t) be the angle system with 

fkis^t)=f^-s, 7f=(s,t)=7;^+s, 

a*^(s, t) = all + "iKS' t) = o^ki - 

and all other angles the same as in q; see Figure 7 (right). Let 

U = {(s, i) G I g(s, t) e A{f, 9, S)}. 

This is a bounded open subset of Let /(s,t) = Ff{q{s,t)). For (so,to) £ U, 
the two tetrahedra (a*i (sq, ^o)) 7*^ (so, ^o)) and (a*^ (sq, io), 7*^ (sO) ^o)) fit together 
iff (soj^o) is a critical point of f{s,t). To see this, apply the same reasoning that 
was used to prove Lemma 1 in Section 4. Also, f{s,t) is strictly concave on U, 
because it is a restriction of the strictly concave function Ff- (Lemma 2) to an 
affine subspace. Finally, the maximum of f{s,t) on the compact set U cannot be 
attained on the boundary dU. To see this, apply the same reasoning that was 
used in Case 1 and in Case 3(a) above. (The tetrahedra (a*^ (s, i), 7*^ (s, t)) and 
(a*^(s,i),7*^(s,t)) cannot degenerate badly for (s,t) e U.) Hence the restriction 
of f{s,t) to U attains its maximum at some (sm,^m) S U, and we have found 
q = q{sm, tm)- Finally we obtain q e A{T, 9, E) with Fr{q) = Ff{q) > Frip) from 
q by flipping the edge kl. 

Since we have thus dealt with the last Case, this completes the proof of the 
Claim, and hence the proof of Lemma 3. 

7. Volume computations. Proof of Lemma 4 

In this section we derive Equation (9) for the volume of a tetrahedron with one ideal 
and three hyperideal vertices, that is, we prove Lemma 4 of Section 4. Vinberg [36] 
[37] derived a formula — in fact, the same formula — for the volume of a tetrahedron 
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D 



a 



C 




'A 



B 



Figure 9. The bircctangular tetrahedron Pi with two ideal vertices A and D. The 
dihedral angles at the unlabeled edges are 

with one ideal and three finite vertices. Here we follow a very similar path. We 
subdivide a tetrahedron with one ideal and three hyperideal vertices into three 
special pyramids (Section 7.5). A volume formula for these is derived in Section 7.4 
with the help of yet other volume formulas that we present in the following sections. 

7.1. The volume of a birectangular tetrahedron with two ideal vertices. 

A tetrahedron with vertices ABCD is called birectangular or an orthoscheme if the 
edge AB is perpendicular on the side BCD and the edge CD is perpendicular on 
the side ABC. Then the dihedral angles at edges AC, BD, and BC are ^. A 
formula for the volume of a birectangular hyperbolic tetrahedron as a function of 
the remaining three dihedral angles was already derived by Lobachevsky [21], see 
also Coxeter [13]. We are only interested in the case of a birectangular tetrahedron 
Pi whose vertices A and D are ideal (see Figure 9). Because the dihedral angles 
sum to TT at the ideal vertices it follows that the angles at AB and CD are equal, 
say, to a, and the angle at AD is ^ — a. Milnor derived the particularly simple 
volume formula 



by direct integration [22] [24]. 

7.2. The volume of an ideal prism. Let P2 be a triangular prism with all vertices 
at infinity. Such a prism is always symmetric with respect to a planar reflection 
that interchanges the triangular faces. (This is so because any three points on 
5'^, the sphere at infinity, can be mapped to any other three points on S*^ by an 
orientation reversing Mobius transformation of S*^, and such a transformation is 
the restriction of a hyperbolic reflection.) Formula (32) below for the volume of P 
is derived in [19]. For the reader's convenience we reproduce the argument. Let 
the interior dihedral angles be a, (3, 7, a', 7', as shown in Figure 10 (left). The 
symmetry plane intersects the side faces of the prism orthogonally in the dashed 
triangle. Hence there exists an ideal prism with dihedral angles a, P, 7, iff there is 
a hyperbolic triangle with these angles, i.e. iff 



Vol(Pi) = -JI(a) 



(31) 



a + /3 + 7 < TT. 



The dihedral angles sum to tt at each ideal vertex, hence 



TT — a — /? + 7 



TT + Q' — /3 — 7 



/3' 



TT — a + /? — 7 



7 



2 



a ~ 



2 



2 
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Figure 10. The ideal prism P2- To calculate its volume, we subdivide it into three 
ideal tetrahedra. 



9 




Figure 11. The tetrahedron P3 with one hyperideal vertex (marked o) and three 
ideal vertices. The points where the edges intersect the sphere at infinity are 
marked •. The tetrahedron is truncated at the polar plane of the hyperideal vertex. 



Lemma 12 (Leibon [19]). The volume of the prism P2 is 
Vol(F2) = JI(a) + JI(/3) + Jl(7) 

+ ji (!L±^_A^) + ji ( "-" + ^"^ ) 

+ Jl C-^'/^^' )+Jl C-"-/-' ). (32) 

Proof. Subdivide the the prism P2 into three ideal tetrahedra as shown in Figure 10 
(middle, right). Since in an ideal tetrahedron the dihedral angles at opposite edges 
are equal [22] [24] , most of the dihedral angles of the three tetrahedra are equal to 
some angle of the prism; see Figure 10 (right). The remaining two dihedral angles, 
A and /U, are obtained by considering how the angles of the prism are sums of angles 
of the tetrahedra: 

A = /3-/3= , ^ = 7r-7. 

Now the volumes of the three tetrahedra are obtained from Equation (24). Take 
the sum and note that JI(7r — a;) = — Jl{x) to obtain Equation (32). □ 



7.3. The truncated volume of a tetrahedron with one hyperideal and 
three ideal vertices. Let P3 be a tetrahedron with one hyperideal and three 
ideal vertices as shown in Figure 11. 
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A 



Figure 12. The special pyramid P4 (left). Decomposition of a tetrahedron with 
one hyperideal and three ideal vertices into one special pyramid and four birectan- 
gular tetrahedra, two of which are mirror symmetric to each other (right). 



Lemma 13. The truncated volume (see Definition on p. 12) of is 



Vol(P3) = ^(jI(a) + Jl(/?) + Jl(7) 



+ Jl(" + "^^"^)+Jl("'" + ^'^) 

+ + (33) 

Proof. If you reflect the truncated tetrahedron at the truncation plane, you get an 
ideal prism. The volume of the truncated tetrahedron is therefore half the volume 
of the ideal prism, Equation (32). □ 

Remark. The same volume formula holds when the apex is finite instead of hyper- 
ideal [36] [37]. In that case a + /3 + 7 > tt. 

7.4. The truncated volume of a special pyramid. Let P4 be a pyramid over 
a quadrilateral base, such that the apex C is at infinity, one lateral edge CD is 
perpendicular to the base, the vertex O of the base that is opposite D is hyperideal, 
and at the other two vertices the angles of the base quadrilateral are ^; see Figure 12 
(left). Let the interior dihedral angles at the edges emanating from O be a, (3, and 
7 as shown. They satisfy a, (3 < tt/2 and a + /? + 7 < tt. Since the sum of dihedral 
angles at the four-valent apex C is 2tt and two of the incident edges have dihedral 
angle ^, the dihedral angle at edge CD is tt — 7. 

Lemma 14. The truncated volume of P4 is 

m) = ^fji(7)+Jif "+":^'^ )+ji^"'"+^'^ 



2V V 2 J \ 2 



/TT — a — /3 + 7\ /TT — a — /3 — 7\, 
-JI +JI . (34) 



Proof. Extend the edges of the base emanating from O until they intersect the infi- 
nite boundary at the ideal points A, B, see Figure 12 (right). The truncated volume 
of the tetrahedron P3 with vertices ABCO is given by Equation (33). It can be 
partitioned into the special pyramid P4, the tetrahedron ABCD, and two birectan- 
gular tetrahedra as shown. The volumes of the birectangular tetrahedra are given 
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C 



1 



D 




O 



Figure 13. The base may also be a self-intersecting quadrilateral. Here, /3 > ^. 



Figure 14. Tetrahedron with one ideal and three hyperideal vertices, partitioned 
into four special pyramids 

by Equation (31); they arc ^ JI(a) and i The tetrahedron ABCD is sym- 

metric with respect to reflection at the plane that contains edge CD and intersects 
edge AB orthogonally. This symmetry plane splits the tetrahedron ABCD into 
two symmetric birectangular tctrahcdra. At edge CD, each of them has a dihedral 
angle of ^(tt — a — /3 + 7). (To see this consider the dihedral angles of all the pieces 
at edge CD; their sum is 2tt. Thus the angle in question is ^{2n — a — (5— (tt — 7)).) 
The volume of the tetrahedron ABCD is therefore JI(i(7r — a — /3 + 7)). Subtract 
from the truncated volume of P4 the volumes of tetrahedron ABCD and the two 
birectangular tetrahedra to obtain Equation (34). □ 

Remark. The same volume formula holds when O is finite instead of hyperideal [36] 
[37] . In that case a + (3 -\- ^ > tt. 

Equation (34) holds also when the base is a self-intersecting quadrilateral; see 
Figure 13. In this case one of the angles a or /3 is greater than |. We have to 
regard the pyramid P4 as a difference of two tetrahedra (with signs as indicated 
in Figure 13), and its volume as the difference of the volumes of these tetrahedra. 
One can derive Equation (34) by a similar construction. 

7.5. The truncated volume of a tetrahedron with one ideal and three 
hyperideal vertices. Finally, consider a tetrahedron with one ideal and three 
hyperideal vertices as shown in Figure 4. Drop the perpendicular from the ideal 
vertex onto the plane of the opposite face and subdivide the tetrahedron into four 
special pyramids (the bases of which may be self-intersecting quadrilaterals) as 
shown Figure 14. The volume formula for the tetrahedron. Equation (9), is obtained 
as the sum of the volumes of the four special pyramids, given by Equation (34). 
(Note Jl(|- — x) = — JI(|- +x).) This concludes the proof of Lemma 4. 
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